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We present measurements of branching fractions and charge asymmetries for seven B-meson 
decays with an n, r\ or to meson in the final state. The data sample corresponds to 89 million 
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BB pairs produced from e + e annihilation at the T(4S) resonance. We measure the following 
branching fractions in units of 10" 6 : B(B + -> ?77r + ) = 5.3±1.0±0.3, B(B + -» r)K + ) = 3.4±0.8±0.2, 
B(B° -> r]K ) = 2.9±1.0±0.2 (< 5.2, 90% C.L.), B(B+ -> r)'-K+) = 2.7±1.2±0.3 (< 4.5, 90% C.L.), 
B(B + -> UJTT+) = 5.5±0.9±0.5, B(B + -> luK+) = 4.8±0.8±0.4, and B(B° -» wK°) = 5.9lij±0.5. 
The charge asymmetries are A c h{B + -> r?7r+) = -0.44 ± 0.18 ± 0.01, A*(B + -> r/K + ) = -0.52 ± 
0.24 ± 0.01, A c h(B + -» LJ7T+) = 0.03 ± 0.16 ± 0.01 and A ch (B+ -» = -0.09 ± 0.17 ± 0.01. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



We report results of measurements of B decays to the 
charmless final states r/K , rjir + , tjK + , rj'ir + , ujK°, ujtt + , 
and loK + Q] . Only the last two of these decays have been 
observed previously [3-[JE| Measurements of the related 
B — > rfK decays were published recently 0. Charmless 
decays with kaons are usually expected to be dominated 
by b — > s loop ("penguin") amplitudes, while b — > u 
tree transitions are typically larger for the decays with 
pions. However the B — > r\K decays are especially inter- 
esting since they are suppressed relative to the abundant 
B — > rj'K decays due to destructive interference between 
two penguin amplitudes |g. Thus the CKM-suppressed 
b — ► u amplitudes may interfere significantly with the 
suppressed penguin amplitudes. This tree-penguin in- 
terference may lead to large direct CP violation in the 
•qK + decay as well as r/7r + , and 7/71" + Q; numerical esti- 
mates have been provided in a few cases ||. We search 
for such direct CP violation by measuring the charge 
asymmetry A c h = (r _ — r + )/(r~ + T + ) in the rates 
= T(B ± — > / ), for each observed charged final state 

f ± - 

Charmless B decays are becoming useful to test the ac- 
curacy of theoretical predictions such as QCD factoriza- 
tion |9j ■ Phenomenological fits to the branching fractions 
and charge asymmetries can be used to understand the 
importance of tree and penguin contributions and may 
even provide sensitivity to the CKM angle 7 (lo| . 

The results presented here are based on data col- 
lected with the BABAR detector at the PEP-II asym- 
metric e + e~ collider 0] located at the Stanford Lin- 
car Accelerator Center. An integrated luminosity of 
81.9 fb _1 , corresponding to 88.9 ± 1.0 million BB pairs, 
was recorded at the T(45) resonance (center-of-mass en- 
ergy y/s = 10.58 GeV). 

Charged particles from the e + e~ interactions are de- 
tected, and their momenta measured, by a combination of 
a vertex tracker (SVT) consisting of five layers of double- 
sided silicon microstrip detectors, and a 40-layer cen- 
tral drift chamber, both operating in the 1.5-T magnetic 
field of a superconducting solenoid. We identify photons 
and electrons using a CsI(Tl) electromagnetic calorimeter 
(EMC). Further charged particle identification (PID) is 
provided by the average energy loss (dE/dx) in the track- 
ing devices and by an internally reflecting ring imaging 
Cherenkov detector (DIRC) covering the central region. 

We select 77, 77', u, K®, and tt° candidates through the 
decays 77 -> 77 (?7 77 ), 77 -> 7r + 7r~7r° (773^), r/' -> rjir+n^ 



Wr,™), V' -> P°l iVp-y), U^TT+TT 7T° , P ° -> 7T+7T , 

Kg — * 7r + 7T _ , and n° — + 77. We make the following re- 
quirements on the invariant mass (in MeV) of their final 
states: 490 < to 77 < 600 for 77 77 , 520 < < 570 for 

773^, 910 < (m^TT-n, m pl ) < 1000 for if , 735 < < 825 

for uj, 510 < m-xir < 1070 for p°, and 120 < m 77 < 150 for 
7T°. For K® candidates we require 488 < < 508, the 
three-dimensional flight distance from the event primary 
vertex to be greater than 2 mm, and the angle between 
flight and momentum vectors, in the plane perpendicular 
to the beam direction, to be less than 40 mrad. 

We make several PID requirements to ensure the iden- 
tity of the pions and kaons. Secondary tracks in 773^, 77', 
and ui candidates must have DIRC, dE/dx, and EMC 
outputs consistent with pions. For the B + decays to an 
77 or u> meson and a charged pion or kaon, the latter (pri- 
mary) track must have an associated DIRC signal with 
a Cherenkov angle within 3.5 standard deviations (cr) of 
the expected value for either a ir or K hypothesis. 

A -B-meson candidate is characterized kincmat- 
ically by the energy-substituted mass ttoes = 
[{\ s + Po ' Pb) 2 /Eq — p 2 B ]z and energy difference 
AE = Eg — \\fs, where the subscripts and B refer to 
the initial T(AS) and to the B candidate, respectively, 
and the asterisk denotes the T(4S) frame. The resolution 
on AE (rn E s) is about 30 MeV (3.0 MeV). We require 
\AE\ < 0.2 GeV and 5.2 < to E s < 5.29 GeV. 

Backgrounds arise primarily from random combina- 
tions in e + e~ — > qq events. We reject these by using 
the angle 9t between the thrust axis of the B candidate 
in the T(45) frame and that of the rest of the charged 
tracks and neutral clusters in the event. The distribution 
of |cos6*t| is sharply peaked near 1.0 for combinations 
drawn from jet-like qq pairs, and nearly uniform for B- 
meson decays. We require |cos#t| < 0.9, for all modes 
except the high-background B + — > r)' pi Tt + decay, where 
we determine that the sensitivity is maximal for a 0.65 
requirement. We also use, in the fit described below, a 
Fisher discriminant T that combines four variables: the 
angles with respect to the beam axis of the B momen- 
tum and B thrust axis (in the T(45) frame), and the 
zcroth and second angular moments Lo,2 of the energy 
flow about the B thrust axis. The moments are defined 
by Lj = J2iPi x l C0S ^| J i where 0j is the angle with re- 
spect to the B thrust axis of track or neutral cluster i, pi 
is its momentum, and the sum excludes the B candidate. 

For the 77 — * 77 modes we use additional event- 
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selection criteria to reduce BB backgrounds from sev- 
eral charmless final states. We reduce background from 
B — > 7t + tt°, K + tt°, and K°ir° by rejecting ry 77 candi- 
dates that share a photon with any ir° candidate hav- 
ing momentum between 1.9 and 3.1 GeV/c in the T(4S) 
frame. Additionally, we require E 1 < 2.4 GeV to sup- 
press background from B — ► K*-f and related radiative- 
penguin decays. From Monte Carlo (MC) simulation 
we estimate that the residual charmless BB background 
is negligible for all decays except those with r\ — ► 77 and 
i]' — > p°7, where we include in the fit described below 
a BB component (which is less than 0.5% of the total 
sample in all cases). 

We obtain yields and A c h from extended unbinned 
maximum-likelihood fits, with input observables AE, 
toes, J 7 , rn ICS (the mass of the 77, 77', or ui candidate), 
for the uj decays, H = |cos0jj|, and for charged modes 
the PID variable S„ t K. The helicity angle Or is defined 
as the angle, measured in the ui rest frame, between the 
normal to the u> decay plane and the flight direction of the 
uj. We incorporate PID information by using S-K 

(Sk), 

the number of standard deviations between the measured 
Cherenkov angle and the expectation for pions (kaons). 

For each event i, hypothesis j (signal, continuum back- 
ground, BB background), and flavor (primary n + or K + ) 
k, we define the probability density function (PDF) 



jk 



x [PM)] [Pjm] ■ 



(1) 



The terms in brackets for S and TL pertain to modes with 
charged track or ui daughters, respectively. The absence 
of correlations among observables in the background Pj k 
is confirmed in the (background-dominated) data sam- 
ples entering the fit. For the signal component, we cor- 
rect for the effect of the neglect of small correlations (see 
below). The likelihood function is 



N 



£ = ex P (-]Ty, fe )n 



j • k 



jk 



j,k 



(2) 



where Yjk is the yield of events of hypothesis j and flavor 
k found by maximizing C, and N is the number of events 
in the sample. 

We determine the PDF parameters from simula- 
tion for the signal and BB background components, 
and from (iues-iAE) sideband data for continuum 
background. We parameterize each of the functions 
^(toes), V sig {AE k ), Vj(f), Vj(S k ) and the peaking 
components oi'Pj(m Tes ) with either a Gaussian, the sum 
of two Gaussians or an asymmetric Gaussian function 
as required to describe the distribution. Slowly varying 
distributions (mass, energy or helicity-angle for combina- 
toric background) are represented by linear or quadratic 
dependencies. The peaking and combinatoric compo- 
nents of the ui mass spectrum each have their own Ti 



shapes. The combinatoric background in toes is de- 
scribed by the function x\J\ — x 1 exp [— £(1 — £ 2 )], with 
x = 2toes/\/* an d parameter £. Large control samples 
of B decays to charmed final states of similar topology 
are used to verify the simulated resolutions in AE and 
toes- Where the control data samples reveal differences 
from MC in mass or energy offset or resolution, we shift 
or scale the resolution used in the likelihood fits. 

In Table [i] we show for each decay mode the measured 
branching fraction, together with the quantities entering 
into its computation. Typically seven parameters of the 
background PDF are free in the fit, along with signal and 
background yields, and for charged modes the signal and 
background A c h- For calculation of branching fractions, 
we assume that the decay rates of the T(AS) to B + B~ 
and B°B° are equal. For the 77 and rj' decays, we combine 
results from the two decay channels by adding the values 
of — 21n£, taking proper account of the correlated and 
uncorrelated systematic errors. The estimated purity is 
the ratio of the signal yield to the effective background 
plus signal; we estimate the effective background by tak- 
ing the square of the uncertainty of the signal yield as the 
sum of effective background plus signal. In Figs. and 
12 we show projections onto toes and AE after requiring 
S„ x ~ 2 [for (a)-(d)] and requiring that the signal like- 
lihood (computed ignoring the PDF associated with the 
variable plotted) exceeds a mode-dependent threshold. 

The statistical error on the signal yield and A c h is 
taken as the change in the central value when the quan- 
tity — 21n£ increases by one unit from its minimum 
value. The significance is taken as the square root of the 
difference between the value of — 21n£ (with systematic 
uncertainties included) for zero signal and the value at its 
minimum. For r]K and rj'ir + we quote a 90% confidence 
level (C.L.) upper limit, taken to be the branching frac- 
tion below which lies 90% of the total of the likelihood 
integral in the positive branching fraction region. For the 
charged modes we also give the charge asymmetry Ach- 

Most of the yield uncertainties arising from lack of 
knowledge of the PDFs have been included in the sta- 
tistical error since most background parameters are free 
in the fit. Varying the signal PDF parameters within 
their estimated uncertainties, we estimate the uncertain- 
ties in the signal PDFs to be 1-3 events. We verify the 
validity of the fit procedure and PDF shapes by demon- 
strating that the likelihood of each fit is consistent with 
the distribution found in simulation. 

Uncertainties in our knowledge of the efficiency, found 
from auxiliary studies, include 0.8iV t %, 2.5iV 7 %, and 3% 
for a K® decay, where N t and 7V 7 are the number of 
signal tracks and photons, respectively. Our estimate of 
the B production systematic error is 1.1%. The neglect 
of correlations among observables in the fit can cause a 
systematic bias; the correction for this bias (< 10% in 
all cases) and assignment of systematic uncertainty (1- 
5%), is determined from simulated samples with varying 
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TABLE I: Signal yield, estimated purity P, detection efficiency e, daughter branching fraction product, significance (including 
systematic uncertainties), measured branching fraction, background (All) and signal (Ach) charge asymmetries for each mode. 
For B° — > r/K° and B + — » t]'tt + , the 90% C.L. upper limit is also given. 



Mode 


Yield 


P (%) 


€ (%) J! &i (%) Signif. 




S(10" 6 ) 


A qq 




??37r7r + 


28^° 


30 


93 


93 


A A 






-0.004 ± 0.010 


-0.52 ±0.31 




59 ± 14 


31 




30 


0.0 




5.2 ± 1.3 


-0.001 ±0.011 


-0.41 ±0.22 


r77T + 














5.3 ± 1.0 ± 0.3 


-0.003 ± 0.008 


-0.44 ± 0.18 ± 0.01 




15l 8 7 


24 


93 


93 






3.i±i:S 


-0.008 ±0.016 


-0.43 ±0.51 


ri 11 K + 


38 ±11 


33 


23 


39 


5.3 




3.5 ±1.1 


-0.011 ±0.016 


-0.55 ±0.26 


T]K + 










6.1 




3.4 ± 0.8 ± 0.2 


-0.010 ±0.011 


-0.52 ± 0.24 ± 0.01 




2.6±|:i 


20 


23 


8 


0.8 




1-8±|;1 








8-6±|| 


47 


24 


14 


3.2 










T)K° 










3.3 


2 


9 ± 1.0 ± 0.2 (< 5.2) 






n' 7T+ 

'IrfTTlT " 




38 


28 


17 


3.9 










77p 7 7T + 






17 


30 






















3.4 


2 


7 ± 1.2 ± 0.3 (< 4.5) 






UJTT + 


101 ± 17 


37 


23 


89 


9.1 




5.5 ± 0.9 ± 0.5 


0.012 ±0.006 


0.03 ± 0.16 ± 0.01 


uiK + 


83 ± 14 


39 


22 


89 


10.0 




4.8 ± 0.8 ± 0.4 


-0.003 ± 0.009 


-0.09 ± 0.17 ± 0.01 




33l 9 8 


51 


20 


31 


7.5 




5.9+11 ± °- 5 








5.2 5.22 5.24 5.26 5.28 -0.2 -0.1 0.1 0.2 



m ES (GeV) AE (GeV) 

FIG. 1: Projections of the B candidate thes and AE for (a, 
b) B + — > 777r + , and (c, d) B + — > r]K + . Points with errors 
represent data, shaded histograms the r\ — > 7r + 7r _ 7r° subset, 
solid curves the full fit functions, and dashed curves the back- 
ground functions (the peaking BB background component is 
negligible). These plots are made with a requirement on the 
likelihood and thus do not show all events in the data samples. 

background populations. Published data [L4J provide the 
uncertainties in the B-daughter product branching frac- 
tions (1%). Selection efficiency uncertainties are 1% (3% 
in B + — > ?7p 7 7r + ) for cos^t and ~1% for PID. Using sev- 
eral large inclusive kaon and B-decay samples, we find 
a systematic uncertainty for A c h of 1.1% due mainly to 
the dependence of reconstruction efficiency on the charge 
of the high momentum charged track. The values of Af! 
(see Table Pi provide confirmation of this estimate. 

In conclusion, we find significant signals for five B- 
meson decays. The measured branching fractions, and 
for the B+ modes the charge asymmetries, are given in 
Table D These are the first charge asymmetry measure- 
ments for the decays B + — > rjir + and B + —* r)K + , since 
these modes along with B° — > ujK have not been ob- 




5.2 5.22 5.24 5.26 5.28 -0.2 -0.1 0.1 0.2 

m ES (GeV) AE (GeV) 



FIG. 2: Projections of the B candidate thes and AE for (a, b) 
B+ -> U7V+- (c, d) B+ -> coK+; and (e, f) B° -» uK°. Points 
with errors represent data, solid curves the full fit functions, 
and dashed curves the background functions. These plots are 
made with a requirement on the likelihood and thus do not 
show all events in the data samples. 



served previously. We quote 90% C.L. upper limits for 
the B° — > r]K° and B + — > rj'n + branching fractions, 
where the significances are only 3.3<r and 3.4<r, respec- 
tively. All branching fraction and charge asymmetry 
measurements are consistent with, but more precise than, 
previous measurements |2|, y, |4J, Il5| . Though uncertain- 
ties are large, the values of A c h for the two decays with u) 
mesons are small as expected theoretically; the consisten- 
cies with zero asymmetry for B + — ► rjTr + (B + — > r]K + ) 
are 2.4<r (2.1a). These are channels in which large asym- 
metries may be anticipated 0. 
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